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Brief Summary 


The principal research effort accomplished under NASA Grant NAG 5-489 
was the research reported in Attachment 1 entitled "Energy Balance Constraints 
on Gravity Wave Induced Eddy Diffusion in the Mesosphere and Lower 
Thermosphere." This paper is in press and should appear in the November 1985 
issue of Journal of Geophysical Research--Atmospheres . 

In this paper we make a strong case for turbulent diffusive transport 
of constituents and potential temperature at rates significantly less than used 
by most workers in the field. Recent analysis of microwave derived water vapor 
mixing ratio profiles in the upper mesosphere by the PI substantiate this 
conclusion of the paper. 

We also argued in this paper that larger diffusion coefficients are 
permitted for chemically active species than the constraints derived from 
thermodynamic considerations. Recent work by the PI shows that the large dif- 
fusion coefficients used by researchers to transport odd oxygen can be supplied 
by gravity waves whose amplitudes are near their saturation value. 

In addition some preliminary research was carried out in non-LTE 
radiative transfer. This research will be continued in subsequent NASA grants 
and used in the analysis of SME data. 
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Abstract 


Breaking gravity waves generate and maintain a background level of tur- 
bulence which is capable of producing substantial cooling and/or heating in the 
upper mesosphere and lower thermosphere. The net thermodynamic effect of 
breaking gravity waves is critically dependent on the eddy Prandtl number (P^) 
applicable to mesospheric turbulence. When P^ -> 1, our calculations of the heat 
budget for the mesopause region imply that the globally averaged eddy or tur- 
bulent diffusion coefficient cannot exceed 10 6 cm 2 s _1 . This upper limit on tur- 
bulent diffusion applies "to^oth potential temperature transport and chemically 
inert tracer transport when radiative damping is negligible. For chemically 
active species larger diffusion coefficients are permitted, because the effec- 

o 

tive eddy diffusion coefficient is increased by an additive term L/2y , where L 
is the chemical loss rate and y is the vertical wavenumber. 

For Pt - 4-6, the turbulent diffusion of momentum (D^) is sufficiently 
greater than the turbulent diffusion of heat (D^) that the conversion of gravity 
wave energy to heat with high efficiency nearly balances the divergence of the 
downward eddy heat flux in the wave breaking zone. Therefore the heat budget of 
the mesopause region would no longer provide a powerful and useful constraint on 
&H* If P t exceeds 6 with high efficiency for energy conversion to heat, gravity 
waves would heat the mesosphere throughout the wave breaking region. 
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1. Introduction 

Lindzen's (1981) parameterization of turbulence and stress generated by 
breaking gravity and tidal waves initiated considerable research on constituent 
transport and the momentum and heat budget of the mesosphere (e.g. Holton 1982, 
1983; Dunkerton, 1982; Schoeberl et al., 1983; Apruzese et al., 1984; Fritts and 
Dunkerton, 1985; Garcia and Solomon, 1985). The identification of a specific 
mechanism which produces turbulent mixing allows a quantitative description of 
associated effects even though it may be highly parameterized. On the basis of 
the research referenced above it would appear that the intensity of the turbu- 
lence as represented by the eddy or turbulent diffusion coefficient, D(cm s ), 
must be greater than 10 6 cm 2 s -1 from studies of the momentum budget and consti- 
tuent transport, whereas our previous study of the heat budget (Apruzese et al., 
1984) concluded that D < 10 6 cm 2 s -1 with a preferred value of 6 x 10 5 cm 2 s -1 . 

Taken at face value we might deduce that the effective eddy Prandtl 
number must be greater than one. (The eddy Prandtl number, P^, is the momentum 
diffusion coefficient, D M , divided by the thermal diffusion coefficient, 

In a turbulent atmosphere tracer transport and potential temperature transport 
should be described by the same coefficient, D^, in the limit of a chemically 
inert species and adiabatic motion, if K theory adequately describes turbulent 
transport processes (see Kraichnan, 1976, for comments on this simplification). 

The studies of Allen et al. (1981) and Garcia and Solomon (1985) clearly suggest 
that Dh exceed 10 6 cm 2 s _l for chemically active species in contrast to the upper 
limit on obtained by Apruzese et al . (1984). In the case of Garcia and Solomon 
(1985) their net had only a weak dependence on the chemical loss rate and thus 
should be approximately equal to the heat diffusion coefficient. They assumed 
Pt = 1 and calculated Dj*j with Holton's (1982) modification of Lindzen's (1981) 
parameterization of breaking gravity waves. As indicated above they obtained 
> 10 6 cm 2 s _1 in the mesopause region, with largest values during solstice. 
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The principal objective of this report is to re-examine the constraints 
on turbulence imposed by the heat budget and determine whether there is suffi- 
cient theoretical evidence to support the hypothesis that the eddy Prandtl 
number is greater than one in the mesosphere. The mesopause thermal structure 
is calculated with turbulent diffusion coefficients commonly used in chemical 
models and deduced from mean zonal wind deceleration. As seen below, extreme 
mesopause temperatures of less than 100 K are produced by the large net cooling. 

2. Model 

The model adopted for this study includes the non-LTE version of the 
Apruzese et al. (1982) two-stream IR algorithm, the physical processes of molecular 
and turbulent heat conduction, odd oxygen transport, gravity wave dissipation, 
and NO IR cooling as described in Apruzese et al. (1984). The approximate IR 
algorithm incorporates exact cooling to space escape probabilities from Fels and 
Schwarzkopf’s (1981) extensive tables on COg transmittances to generate 
accurate cooling to space rates. Radiative exchange between atmospheric layers 
is treated approximately; consequently only thermal disturbances with vertical 
wavelengths greater than 8 km can be treated with sufficient accuracy in the 
model. None of the results presented below are sufficiently sensitive to our 
approximate IR cooling and heating rates to warrant a substantially more 
accurate treatment of IR radiative transfer at this time. 

Vertically propagating gravity waves carry' energy which is preferen- 
tially deposited as heat In the wave breaking zone as a result of the turbulence 
created. Schoeberl et al. (1983) gave a simple analytic expression for this 
heating in terms of the efficiency e with which gravity wave energy is converted 
to heat 

N 2 D 

“,-iri 1 * p t K 5 

P 


( 1 ) 



where N is the buoyancy frequency, and Cp is the specific heat at constant 
pressure. Also turbulence generated by breaking gravity waves leads to a down- 
ward turbulent heat flux whose divergence (convergence) results in cooling 
(heating). The usual expression for this process acting on the mean temperature 
field is 



where > 0 for heating, p is mass density, T is temperature, g is gravita- 
tional acceleration, and z is altitude. 

Given the variety of physical and radiative processes outlined above 
and discussed in depth in Apruzese et al . (1984) the globally averaged 
equilibrium temperatures in the mesosphere and lower thermosphere were calcu- 
lated as a function of the eddy diffusion coefficient, D^, and eddy Prandtl 
number, P*. In computing these temperatures it was presumed that the correla- 
tion between the turbulent diffusion and temperatures could be represented by 
the product of their global averages. In other words if <> represents a global 
average then we assumed 

< _i D _3 T > _J_ <D> _3 < T > 

3Z 3Z dz 3Z 

We have tested this assumption under solstitial conditions using an analytic 
model of the mesospheric wind and temperature structure. Using Lindzen's for- 
mula D « (IT-c) 4 we found that using global averages to represent the heat flux 
is accurate to better than 1%. 
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The global heat budget is also affected by the convergence of the down- 
ward heat flux associated with the anticorrelation of the zonal mean vertical 

velocity, w, (a result of gravity wave breaking) and the zonally averaged tem- 

1 a 

perature, i.e. the term ~ < p wT> in the thermodynamic energy equation. This 
term is substantially smaller than the net heating term due to gravity waves, H g 
+ Cg. For example at the mesopause under solstitial conditions the model of 
Garcia and Solomon (1985) with a mean meridional circulation driven by gravity 
wave breaking predicts < D > = 1.5 x 10 6 cm s -1 and < wT > = 16 K cm s 1 . The w 

I/O 

and T fields have an approximate height dependence of p . The convergence of 
this "1reat flux is spread over at least 2 or 3 scale heights and has a magnitude 
of at most 10% of Hg + Cg. This is not surprising as w and T have large ampli- 
tudes only at polar latitudes; the global average is accordingly small. 

3. Results 

In our initial set of calculations, four representative turbulent dif- 
fusion coefficients were constructed and are shown in the top panel of Fig. 1. 
Model A is our "best fit" profile from Apruzese et al . (1984) (their model 4), 
although no real effort was expended to generate a truly excellent fit. Model A 
has an asymptotic value for D of 6xl0 5 cm 2 s -1 at high altitudes and is charac- 
terized by very low values throughout the mesosphere. 

Model B is an average of the summer and winter diffusion coefficients 
estimated by Lindzen (1981) for breaking gravity waves. Model C includes, in 
addition to gravity wave generated turbulent diffusion, Lindzen's estimate of the 
contribution from the breaking diurnal tide. These two profiles should be 
regarded as illustrative rather than rigorous globally averaged values of the 
turbulent diffusion coefficient. Finally Model D is the envelope of the highly 
structured vertical eddy diffusion coefficient used by Allen et al . (1981) in a 
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ID model of mesospheric chemistry. Most ID models of mesospheric chemistry 
require eddy diffusion coefficients in excess of 10 cm s and possibly as 
high as 10 7 cm 2 s _1 . With a peak value of 2xl0 6 cm 2 s -1 , Model D is certainly not 
an extreme case. 

In the middle panel of Fig. 1 the four globally averaged temperature 
profiles calculated with their respective turbulent diffusion coefficients are 
compared with the CIRA 1972 temperature profile, which is representati ve of the 
observed globally averaged temperature. In these calculations it was assumed 
that Pt = 1, i.e. D^ = D^ = D, and e = 1. The assumption of a high efficiency 
for conversion of gravity wave energy to heat^is based on numerical simulations 
of Fritts and Dunkerton (1984). Note that Model A underestimates the temperature 
by up to 20°K between 65 and 75 km, the region where IR cooling is dominated by 
GOfehot bands and isotppic bands (Dickinson, 1973; Wehrbein and Leovy, 1982). A 
comparison of our cooling rates with the very recent computations of Dickinson 
(1985) indicates that our rates are too large in the 63-82 km region. As will be 
shown below a 40% reduction in the C0 2 cooling rate eliminates the discrepancy 
between observed and calculated temperatures. The departure of the Model A 
calculated temperature profile above 95 km from CIRA 1972 is due to the uncer- 
tainty in the C0 2 (v 2 ) deactivation rate by 0 collisions as discussed at length in 
Apruzese et al. (1984). The thermosphere is not of fundamental importance in 
this study and is not discussed further. 

In the lower panel of Fig. 1 the divergence of the turbulent heat flux 
is illustrated. For Model A below the mesopause the cooling rate does not 
exceed 5 K d -1 . With Models B and C which have a peak value of D -» 3xl0 6 cm 2-1 , 
the respective mesopause temperatures are 65 and 46 K. Note that at the meso- 
pause the value of D in Model C is significantly larger than Model B’s value. 
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The reason for such cold mesopause temperatures is the extreme divergence of the 
turbulent heat fluxes. In Model B it peaks at - 30 K d _1 , whereas in Model C it 
approaches 350°K d -1 . With solar heating rates of less 5 K d -1 , the equilibrium 
temperature must become very cold to minimize IR cooling and maximize con- 
vergence of downward heat flux from the thermosphere carried by molecular con- 
duction. Model D also yields a cold mesopause (•'100 K) but up to 65 km is in 
close agreement with Model A, which suggests that below 65 km the atmospheric 
heat budget is insensitive to values of D that do not exceed the Model D pro- 
file. Above 80 km the large values of D in Model D generate a divergence of the 
turbulent heat flux in excess of - 63Kd«L* 

To analytically illustrate the strong sensitivity of temperature to 


turbulent diffusion, we write the net heating due to gravity waves as 


H + C = 
9 9 


% 1 + P ^( M ] KS ' 1 


(3) 


when the small term proportional to — y is neglected, Hp 1= ancl C P 

= g(*L + = 

T[zz Cp ) 

3.5 x 10~ H s”*, H » 6 km, Dm in units of 10°cm*s 1 , with P^ = 1 and e = 1 as 


for a diatomic gas. For typical mesospheric values N t ^34 

■\— *♦_— 2 1 1 c 1 n j. * j. _ 1 r»k 2 . — 1 


assumed in Fig. 1 


H g + C g = 7.5 ^1.4 - 1^ D M (10 6 cm 2 s _1 ) K d -1 


(4) 


Thus for constant Dm, Hg + Cg * -7.5 K d" 1 when Dm = 10 6 cn^s -1 or -22 K d _1 for 
Dm = 3xl0 6 cm 2 s _1 , which may be compared against a maximum solar heating rate in 
this region of 5 K d" 1 . Thus when H g + C g < -5 K d -1 (Dm > 7xl0 5 cm 2 s _1 ), the meso- 
pause must cool down to decrease IR cooling and increase convergence of the 

H 

downward thermal conduction heat flux. Note that when xr - 0.7, i.e. Dm « 

••u 

exp(z/1.4H), then Hg + Cg 0. If Dm increases more rapidly with altitude than 
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this, gravity waves will produce a net heating of the region. Similarly if 
increases more slowly with height than exp (z/1.4H) net gravity wave cooling 
occurs. 


Returning to Fig. 1, it is apparent why Model A with an asymptotic 
value of Dm = 6x10 s cm 2 s _1 is close to the upper limit on D^, beyond which 
disastrous consequences on the mesopause temperature result from gravity wave 
cooling. Model D has increasing as exp (z/2H) from 50 to 80 km and hence 
generates net gravity wave cooling over this region. Only when > 10 6 cn^s -1 , 
does this cooling rate become excessive and the mesospheric temperature plunges. 
In Models B and C, D^ > 10 6 cm 2 s -1 from 50 to 75 km and does not increase 
with height more rapidly than exp(z/2H). The large gravity wave cooling rate in 
this region drives the mesospheric temperatures below 100°K above 75 km. 

Examination of Eq. (3) reveals that gravity wave cooling is dominated 
by the last ternfp-^. Although over limited regions Hq can be less than 
1.4H, eventually must level off (Hq » H) and the last term overwhelms the 
other terms with net gravity wave cooling. Only in the limit of large P^ can 
the contribution of this last term be substantially reduced over the entire 
mesosphere. For constant Dm(Hq = ®) when P^ =6, e = 1, then Hg + Cg - 0. With 
finite Hq, a lower Prandtl number will yield negligible gravity wave cooling. 

On the other hand, for large P^(>10) gravity wave breaking will only produce net 
heating. 


To test the effect of the Prandtl number dependence numerically, some 
additional calculations were performed. In these calculations the C0 2 cooling 
rates between 63 and 82 km were reduced by 40% to bring them into closer 
agreement with Dickinson's (1985) more accurate computations and to demonstrate 
that our conclusions do not depend critically on the IR cooling rates. In addi- 
tion the C0 2 (v 2 ) deactivation rate by 0 collisions was adjusted to 
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2x10” 1 1 exp ( -800/T ) cm 3 s _1 (slightly larger than the value adopted by Gordiets et 
al. (1982)) to bring the calculated temperatures in the thermosphere into 
agreement with C1RA 1972. For the reader interested in an extended discussion 
of the issues raised by this adjustment. Section 3.6 of Apruzese et al . (1984) 
should be consulted. 

In Fig. 2 Model 1 contains the above modifications with e = 1, Pt = 4 
and a Dm profile that is - 3.3 times the Model A diffusion coefficient in 
Fig. 1. (The precise expression for Dm is Dm = D 0 D ot> /(D 0 +D„), where D 0 = 3 x 10 5 
exp [(z-65 km/H] cm 2 s _1 , H = scale height, and D*, is the asymptotic value at 
high altitudes given in figure caption.) The asymptotic value of is 5xl0 5 
cm 2 s -1 for Model 1, slightly less than for Model A. The calculated mesospheric 
temperatures in the 65 to 75 km region are now in good agreement with CIRA 1972. 
Also the lower thermospheric temperatures above 95 km are also in good agreement 
with CIRA 1972. The principal region of disagreement is between 80 and 95 km 
where our calculated temperatures exceed CIRA 1972 by up to 15° K, and would 
suggest that enhanced gravity wave cooling can be accommodated by lowering P^. 

If P t = 1, then the Model 2 temperature profile is obtained with a mesopause 
temperature of 140 K. Comparison of Models 1 and 2 illustrates the extreme sen- 
sitivity to Pt that we would expect from our discussion of Eq.(3). It is con- 
fined to the region above 80 km where gravity wave cooling exceeds solar heating. 
In Model 3, Pt * 4 but Dm is twice its value in Model 1. Above 100 km, the temper 
ature is only lowered by 10 K with this factor of 2 increase in Dm* Below 100 
km the temperature change is insignificant. The weak dependence on Dm (or D^) 
when Pt = 4 is due to the large cancellation of gravity wave heating and cooling 
terms. For example at 90 km, Hg + Cg = -2.7 K d 1 for Model 1 and -3.4 K d 1 
for Model 3, with Cg = -6.3 and -10.6 K d -1 , respectively. When compared with 



9 


the solar heating rate of ' 11 K d” 1 , which is rapidly increasing with height at 
90 km, the 0.7 K d” 1 change in Hg + Cg is insignificant. 

4. Discussion and Summary 

The present study clearly demonstrates the importance of the Prandtl 
number for mesospheric turbulence. For large Pt('6) and e » 1, gravity wave 
heating approximately balances gravity wave cooling and the globally averaged 
atmospheric temperature profile is controlled by radiative processes. For P^ > 

6, gravity wave breaking produces only net heating. In the other limit P^ 4 1, 
gravity wave breaking is generally characterized by very strong cooling but with 
small heating at the bottom of the breaking zone. In this limit, the cooling is 
sufficiently strong to impose a powerful constraint on the maximum permissible 
level of turbulence in the mesosphere without plunging the mesopause temperature 
down to 100 K. With allowance for uncertainty in our radiative transfer model 
this upper limit is D = 10 6 cm 2 s -1 . Our preferred asymptotic value is (5-6 )xl0 5 
cm 2 s _1 when P t = 1. Below 70 km, the D profile from Model D is an absolute upper 
1 i mi t . 

From a close examination of Figs. 5 and 6 in Apruzese- et al. (1984), 
one concludes that deducing a local upper limit to D independently of the total 
D profile is risky below the mesopause. For example. Model 3 in Apruzese et al. 
(1984) almost satisfies the upper limit constraint <10 cm s (D^ = 

1.2x10* cm 2 s -1 for high altitudes) and is less than the diffusion profile of 
Model D below 70 km. But the calculated mesopause temperature is - 130 K. 

As discussed above, more rapid transport of constituents is needed than 
allowed by the heat transport. Schoeberl et al. (1983) extended the parameteri- 
zation of linear wave chemical transport by eddy diffusion given in Strobel 
(1981) to include the effect of breaking gravity waves. Their Eqs. (27) and 
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( 28 ) may be combined to yield an effective vertical eddy diffusion coeffcient of 

Kzz '= °H + + | d H 

where y is the vertical wave number and L is the chemical loss rate is units of 
s“*. The first term represents diffusion by the background turbulence created 
by the breaking gravity wave and the second term represents eddy transport by 
the gravity wave for chemically active (L*0) species. Chao and Schoeberl (1984) 
have demonstrated that the last term, which arises from direct diffusive 
transport by the breaking gravity wave, should be negligible rather than 1/2 
because the perturbation vertical velocity and mixing ratio are~almost in 
quadrature. The eddy transport of heat and constituents by a gravity wave in 
the breaking region cannot be represented by a linear diffusi^tTroces’s 
according to their analysis. Schoeberl et al . (1983) showed for chemically 
active species In the mesosphere and lower thermosphere that the second term in 
Eq.(5) can equal or exceed the first term when obeys our upper limit 

constraint. The authors are not aware of any ID chemical model that has tested 

— 62-1 - 
whether K 1Z with <10 cm s can adequately account for constituent transport 

in the mesopause region. If it can be proven that Pt = 1 for mesospheric tur- 

bulence, then the L/2y z contribution to *zz with Dh < 10 cm s would be essen- 

tial. Note that expression (5) is valid only when k(u-c) » y D^, and L, where 

k Is the horizontal wave number, u is the mean zonal wind, and c is the zonal 

phase speed. 

The determination of the eddy Prandtl number in the mesosphere is an 
exceedingly difficult experimental measurement. To date only inferences are 
available that are not definitive or compelling. Although Justus (1967) deter- 
mined P t to be '3 from photographic tracking of rocket released chemical clouds 
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and analysis of turbulent wind data in the 90 - 110 km region, there was con- 
siderable scatter in his data points. Recently Fritts and Dunkerton (1985) have 
theoretically examined constituent .and heat fluxes driven by localized gravity 
wave breaking where the breaking zones are small in vertical extent in com- 
parison to the vertical wavelength. In their analysis this localization of tur- 
bulence leads to an eddy Prandtl number given by 

p t 1 = 1 -|ff + r (6) 

where a is a measure of supersaturation of the breaking gravity wave (a = 1 for 
saturation), n is a measure of the localization of turbulence, i.e. 

D = and ^ = kx+mz-ct is the wave phase with the usual definitions for 

these quantities. They recommend a = 1.3; thus for n = 1.5, = 3.5 and n = 

2, P t - 7. With a = 1, n =1.5 and 2, the respective values of Pt are 3.3 and 
6. With even modest localization of turbulence (n = 1), Pt - 2 for a - 1-1.5. 
The analysis of Fritts and Dunkerton (1985) strongly suggests that Pt * 2 unless 
large supersaturation (a ^ 2} results in the mesosphere from breaking gravity 
waves. Then the eddy Prandtl number can be reduced to s 1. Observations 
by Theon et al. (1967) and Philbrick et al . (1983) indicate supersaturation by 
up to a ' 1.3, but not 2. 

The Fritts and Dunkerton (1985) model is, however, an analytic treat- 
ment of convection effects for a single breaking gravity wave discussed by Chao 
and Schoeberl (1984) and may not accurately represent the breakdown of a gravity 
wave packet. For example, the turbulent region, initially set up by the folding 
of the potential temperature surface, will expand in time as stable air is 
entrained into the turbulent zone. Thus the heat flux will vary in time and a 


12 


turbulent zone which best represents the time mean heat flux may be quite dif- 
ferent from the zone where breaking initially begins. The actual question of 
heat transport by breaking gravity waves may not be settled until detailed 
numerical simulations are performed. 

Finally, it is important to recall that the arguments put forth by Chao 
and Schoeberl (1984) and Fritts and Dunkerton (1985) for heat transport induced 
by a single breaking wave do not include heat transport produced by stable waves 
encountering a turbulent zone and/or being radiatively damped. The effect of 
these processes gives rise to an eddy Prandtl number of only 1. This discussion 
raises the possibility that Pt may indeed be variable with an initial value that 
is large at the onset of wave breaking and as the turbulence develops dissipa- 
tion of stable waves drives the eddy Prandtl number down to an asymptotic value 
Of 1. 

Although this paper has been concerned principally with the globally 

averaged eddy diffusion coefficient, the constraints on it discussed above are 

also applicable locally. For example, the maximum values of occur above the 

cores of the zonal jets near the stratopause. If P^ were - 1,. then a secondary 

circulation would have to be set up in the atmosphere to locally balance the 

6 2—1 

strong cooling when D^ > 10 cm s . 
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Figure Captions 

« 

Fig. 1 (Upper panel) Eddy or turbulent diffusion coefficients described in 
the text as a function of altitude and used to calculate temperature 
profiles (middle panel) which are compared with CIRA 1972 temperature 
profile. (Lower panel) the cooling or heating rate associated with 
the divergence or convergence of the downward eddy heat flux. All 
calculations with P t = 1. 

Fig. 2 Calculated temperature profiles with modified C0 2 cooling rates as a 
function of altitude compared with CIRA 1972 temperature profile which 
illustrate sensitivity to P^. The respective values of P^ and asymp- 
totic values of D^j (D*,) are: (Model 1) [4, 2 x 10 6 ], (Model 2) 

[1, 2 x 10 6 ], and (Model 3) [4, 4 x 10 6 cm 2 s -1 ]. 
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